home *** CD-ROM | disk | FTP | other *** search
/ SGI Freeware 1999 August / SGI Freeware 1999 August.iso / dist / fw_plotutils.idb / usr / freeware / share / ode / limitcycle.ode.z / limitcycle.ode
Encoding:
Text File  |  1998-10-28  |  722 b   |  24 lines
  1. # This example shows a nonlinear system with a stable limit cycle.
  2. # You may run it by typing
  3.  
  4. # ode -f limitcycle.ode | graph -T X -C -x -2 2 -y -2 2
  5. # step 0,2
  6. # step 2,4
  7. # step 4,6
  8. # step 6,8
  9. # .
  10.  
  11. # An orbit converging on the limit cycle will be traced out incrementally.
  12. # If you are using a color X Window System display, each segment of the
  13. # orbit will be a different color.  This is a feature provided by `graph',
  14. # which normally changes the linemode after each dataset it reads. If you
  15. # do not like this feature, you may turn it off by using `graph -B' instead
  16. # of `graph'.
  17.  
  18. yone' = yone + ytwo - yone*(yone*yone+ytwo*ytwo)
  19. ytwo' = -yone + ytwo - ytwo*(yone*yone+ytwo*ytwo)
  20.  
  21. yone = 1
  22. ytwo = 2
  23. print yone, ytwo
  24.